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MAGNETOHYDRODYNAMIC SHOCK STRUCTURE IN A MOVIN G 
MAGNETIC FIELD 

ABSTRACT 

The shock structure around a cylinder immersed in a cold collisionless 

plasma stream containing a weak magnetic field has been investigated in two 

dimensions. The stream pressure on the boundary of the sheath is proportional 
2 

to cos \Jr where f is the angle between the stream velocity and the normal to 
the sheath. The thickness of the sheath and particle density in the sheath 
on the upstream side of the cylinder have been found to be nearly constant. 

A tear drop shaped vacuum is formed downstream from the cylinder. The plasma 
electrons have essentially all the ion stream kinetic energy in the transition 
sheath upstream of the cylinder but the ion stream energy is regained down- 
stream from the cylinder. The importance of this shock structure when applied 
to the earth’s magnetosphere and the generation of Type I comet tails is 
discussed briefly. 

I. Introduction 

The solar wind, consisting of a cold plasma flowing constantly at high 

1,2 

velocity away from the sun, is now known to contain a weak magnetic field 
oriented somewhat perpendicular to the stream velocity vector. Though the 
magnetic field carried along in and moving with the stream velocity is so weak 
that the magnetic field pressure is only one percent of the stream pressure, the 
field prevents the particles from moving independently of one another. The 

result is that any object larger than a few ion gyroradii will create a 

3-<f 

magnetohydrodynamic shock when immersed in the solar wind. 
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The distinctive features of shocks generated in the solar wind are caused 
by the shape of the object creating the shock and the effect of arresting the 
velocity of the magnetic field lines which initially must move with the wind 
velocity since the stream pressure is overwhelmingly greater than the magnetic 
pressure. The plasma and its field are compressed on the upstream surface of 
the obstacle and much of the stream kinetic energy is diverted to gyromotion 
perpendicular to the compressed magnetic field. A simple extension of ordinary 
M theory of the magnetic pinch effect shows that on entering the compressed 
field region through the charge separation electric field on the surface, the 
plasma electrons will acquire essentially all of the kinetic energy of the ion 
stream in the solar wind. This energy is returned to the ion stream after the 
plasma in the compressed field region moves around to the downstream side of 
the obstacle. 

In the work presented here a simple model is discussed in which a cylinder 
of intense magnetic field parallel to its axis presents an obstacle to a cold 
plasma stream containing a weak magnetic field parallel to the cylinder axis. 
Besides the inherent interest of the model in understanding magnetohydrodynamic 
shock structure under such unusual conditions, the model illuminates some gross 
aspects to be expected in astrophysical problems which occur, for example, in 
the boundary of the earth's magnetic field and in the tail rays of comets. 

II. Physical Characteristics of the Transition Sheath Surrounding a Cylinder 
in a Magnetized Plasma Stream 

The simple idealized model considered here is illustrated in Fig. 1. 

and 

A proton-electron plasma, collisionless /with no thermal motion, is streaming 
parallel to the polar axis with velocity v , and encounters a cylinder of 
radius R containing an arbitrarily intense magnetic field sufficient to exclude 
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the incident plasma stream, 
spatially varying magnetic 


(In general the obstacle need jonly^ ^furnish^a 


t 

field pressure, but a cylindrical shape 


has been introduced at the start to simplify the discussion of the paramount 

physical features resulting from the analysis.) Finally, the incident plasma 

parallel to the cylinder axis, 

stream is assumed to contain a weak magnetic field^/which is carried along with 

the plasma. The stream magnetic field exerts a pressure of only one percent 

of the stream pressure and the field traveling with the stream velocity ^v Q is 

maintained by the small gyromotion of the plasma in the moving velocity frame 

of the plasma stream. In this connection it might be noted that the preliminary 

analysis by Dr. M. Neugebauer of the Mariner R plasma data yielded an effective 

plasma thermal energy of the order of 10-30 electron volts in a stream energy 

Dr. Neugebauer^ 

of the order of a kilovolt. /analysis was made by examining the energy 

spread parallel to the stream motion, but the presence of a magnetic field and 
the gyromotion it causes guarantees that the result is valid for the particle 
energy component perpendicular to the stream motion as well. 

The particle pressure on the outer boundary of the transition sheath is 
given by 


2 — „ 2 


Pp = ( n 0 V 0 C0S 1 f ) ( m p V 0 C0S ^) = T1 0 ri ^ r o COS ^ 


( 1 ) 


where n^ is the proton density in the stream m^ is the proton mass and f is 
the angle between the nor ual to the outer boundary surface of the transition 
sheath and the stream velocity vector. The second parenthetical factor 
represents the momentum change of the incident particles at the boundary as they 
enter the sheath. The external pressure confining the particle and enhanced 
magnetic field pressure in the sheath acts perpendicularly to the surface. 

Hence there is no change in particle center of motion momentum component parallel 
to the boundary which would result in a momentum change factor not linearly 
dependent on cos i]r. 
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There is, of course, a very important change in the particle motion 
on crossing the boundary of the sheath and that is that near the stagnation 
point, in particular, much of the stream motion is diverted to gyromotion 
perpendicular to the magnetic field. Although the magnetic field is too 
weak to contribute effectively to the pressure near the stagnation point, 
it does cause the charges to rotate about the lines of force. Particle 
pressure, however, dominates every other consideration and the guiding center 
of the particles are caused to move towards decreasing pressure gradients 
"dragging" the magnetic field lines along with them. The particle pressure 
within the transition sheath is given by 


1 2 1 ,2 2 , 
p = — n m v = — n m (v - v ) 
r p 2spp 2 s p x o m 7 


( 2 ) 


where n g is the ion density within the sheath, v^ is the particle velocity 
component perpendicular to the field, and is the velocity of the guiding 
center, that is, the velocity of the magnetic lines of force. The factor 1/2 
arises because of the two dimensional motion of the particles within the sheath 
as opposed to the monodirectional motion of the stream pressure. In three 
dimensions , where the particle motion is also along the lines of force, the 
factor might be 1/3. 

We thus obtain an expression for the pressure equality on the boundaries 
of the sheath. 


n m v ^cos^ijf + h2/8ji = ~ n m (v ^ - v + (n ^/n ^)H^/8rt = H ^/8n (3) 

opo T o 2sp v o m ' ' s o ' o c v ' 

where H c is the magnetic field intensity in the cylinder at the inner boundary. 

Since the plasma is collisionless , the magnetic field intensity is directly 

proportional to the particle density, and thus the magnetic field intensity 
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anywhere in the sheath is (n /n )H . Particle thermal pressure would be easy 

S u ° the square of the thermal velocity^ 

to include for similar reasons. Since n is proportional to H and/v^, is 

proportional to H if the motion is adiabatic, the particle thermal pressure 

is a constant multiple of the magnetic pressure and may be included by 
2 

multiplying H by an appropriate factor. 

P, 

The change in the external pressure, 4s a function of the distance along 
the boundary of the sheath, 1 , results in a pressure gradient which accelerates 
the plasma within the sheath. 


dv 


m 

n ra = n m v 

s p dt s p m 


dv 
m 

dl 


di + 


d_ 
d 1 


H 

c 

8it 


m 

+ -£. 
1 n 


/ 2 ; 
(v -v 
o m 


, dn 

\ s 

) TT 


(4) 


o 

which may also be obtained by differentiating Eq. 3 with respect to 1. 

2 

Throughout most of the upstream side of the cylinder H is negligible compared 
with the particle pressure and Eq. 3 may be approximated by 


n = 2P/m (v 2 - v 2 ) 
s p o m 

which satisfies Eq. 4 since the last two terms in Eq. 4 are negligible in 

these circumstances c Thus, v o sin\|r throughout the transition sheath on 

the upstream side and n is nearly constant and equal to 2n . As \jr approaches 

s o 

rt/2, however, the magnetic pressure becomes significant and n decreases with 

s 

1 depending on as 


11 2 2 2 
n I 1 + 7- (8fln m v /H ) cos if 

s ou 4 11 * * V opo o 7 Y 

- (3/32)(8nn o m p v o 2 /H o 2 ) 2 cos 4 ^lJ 


(5) 


Using the approximations valid on the upstream side of the cylinder the 


position of the outer boundary to the transition sheath is easily obtained 
by requiring particle conservation. The number of particles arriving per 
second per unit cylinder lengti within a distance, y^, from the polar axis is 
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n Q y 2 v o ; the number leaving the transition region per second per unit cylinder 

length is n (r_-r,)v sin \|r (see Fig. 2). 
s 1 1 o 

r 2 “ r l / ’^ r 2 sin 0 2^ 2 Sin ^ ( 6 ) 

For a cylindrical obstacle for which r^ is constant, © 2 = , and 

r 2 " 2r i (?) 

If the obstacle were not cylindrical a differential form of the constant 



which yields Eq. 7 in the region where n = 2n if r 1 and r„ are now regarded 

sol l 

as functions of 8. If n is a decreasing function of 0 ? as it is near 0 

s 

i|f A*nl2, r 2~ r l i ncreases faster than r 2 and is less than unity. In this 
region, the parametric equations 5 (or rather 3) and 8 must be solved 
numerically. 

The magnetic pressure cannot be neglected in the region 



Fig. 3. 
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III. The Shock Produced Energetic Electrons in the Transition Sheath 

The preceding section has given us a general description of our simple 

two dimensional model problem as illustrated in Figs. 1-3. A weakly 

magnetized plasma moving with velocity v enters a transition sheath where 

the velocity of the particle's guiding centers, that is, the velocity of the 

magnetic lines of force, changes abruptly at the outer boundary of the 

10 11 

transition sheath to a lower velocity ^v^. Ferraro and especially Dungey and 
Rosenbluth have shown how a charge separation layer will be created at any 
boundary where the magnetic field is suddenly increased on the side away 
from the plasma. Beard^* Grad,^ and Sestero^have further 

developed the theory of the electric field on the boundary which serves to 
energize the electrons at the expense of the positive ions. The electric 
field arises due to the difference in momentum between the electrons and 
the positive ions, and thus there is no field in any region where the momenta 
of the opposite charges are equal. 

In the simple two-dimensional problem discussed in this paper, it might 

be anticipated that when the plasma enters and passes through the transition 

sheath the charges will have gained or lost potential energy (depending on 

the sign of their charge) as though they had passed through permeable 

condenser plates at the boundary. Hence, if there is no reverse electric 

field, the electrons will have acquired essentially all of the initial stream 

kinetic energy of the positive ions. 

12 13 

Following Rosenbluth f s aiia lysis in which the plasma did not enter the 
region of increased field intensity, we write, for the changed conditions 
of this problem, the equations of motion of a charge initially moving parallel 
to the x axis 


= ^ E( x) -f -3— Hy 

m v / me 

= - H( x- v ) 

me v nr 


i # 

x 


M 

y 


(11) 

( 12 ) 
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where the coordinate system is illustrated in Fig. 4, q and m are the charge 
and mass of the charged particle, E is the electric field perpendicular to 
the boundary which is tangent to the y,z plane, H is the magnetic field and is 
parallel to the Z axis, and is the velocity of the guiding centers (i c e. 
the magnetic lines of force) which are assumed to be moving entirely parallel 
to the x axis for this part of the analysis. Eq. 12 is easily integrated and 
yields the y component of the particle velocities. 


y = 


me 


ft ( x- 
0-00 


v ) dt 1 
m 


(13) 


Hence the y component of the particle momenta is the same for both electrons 
and protons. 

The particle densities are inversely proportional to the x component of the 

average velocity of the particles. To avoid the impossibly large electric 

fields resulting from excessive charge separation the x component of the 

average velocities of the particles must be kept equal everywhere. This 

requires that in being slowed from an initial velocity v q to the final guiding 

center velocity of V 3 the x-component of the deceleration of the electrons 

m 

must be equal to the deceleration of the ions. Thus Eq. 11 and 13 yield the 
equality 


t *<*> - <^) 2h f 


4 


H (x -v m )dt' = - e_ E( x ) - ( m e 
m. 

—CO 1 




H ( x“V )dt' 


m' 


where the subscripts e and i refer to electron and positive ion and = -ejq.^ = e. 
Neglecting the right hand side of this equation as being of order nWm^ smaller 

than the left hand side we obtain the familiar M theory result 

ft 


(x 
e 


m 


( 14 ) 
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If Eq. 13 and 14 are inserted into Eq. 11, then the integral of Eq. 11 
after multiplication by x is to zero order in m^/nu 


-f 


* 2 2 2e , , 

x, = v - — I E( x )dx' 


Let represent the position at which the incident proton first reaches 

the final guiding center velocity, V , (i.e. the position beyond which the 
compressed field is constant .) then. 


1 “i 


(v/ - O = 


m 



( x ')dx' 


(15) 


Eq. 15 has the very obvious physical interpretation that the total 
potential drop for -oo < x < x is equal to the change in the ion stream 
kinetic energy divided by the ion charge. For x > charge neutrality does 

• f 

not require x^ m only their averages, V^, must be equal. 

Since the momenta of the electrons and protons are equal, no further charge 

separation takes place; the gyroradii of the particles are equal. The momentum 

1/2 

of an ion in the stream m.v has been reduced to (m /m.) m.v in exerting 

x o e l x o 

the pressure on the boundary needed to confine the enhanced magnetic field 

interior to the boundary. The kinetic energy of the ions has been reduced, 

1 2 1 2 1 2 

however, from -r m.v to — m v + — m.v.. . The considerations of the previous 
5 2io 2 e o 2 l M 

section are still valid since the total particle pressure in the transition 

layer is not affected by the electrical potential at the boundary since the 

1 2 1 2 

electron energy in the transition layer is — m.v - m v • 

2 x o 2 i M 

A free hand drawing of the particle motion as it crosses the boundary 
is shown in Fig. 5 for conditions near the stagnation point. The conclusions 
are essentially urchanged at other points on the boundary as illustrated in 
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Fig* 6 if one considers only the component of the guiding center motion 
perpendicular to the boundary. The component of the stream velocity parallel 
to the boundary is the same on both sides of the boundary and therefore the 
problem reduces to the stagnation point case in a coordinate frame which moves 
with the velocity V Q sin \|r parallel to the boundary and v^cos^jr supplants v q in 
the stagnation point considerations. In accelerating to nearly the free stream 
velocity again in its motion away from the stagnation point, the motion of the 
plasma in the transition sheath causes the initial ion kinetic energy to be 
recovered at the expense of the electron energy. 


Cone lusions 

A cold plasma stream containing a weak magnetic field (field pressure much 

less than particle stream pressure) will be compressed on the upstream surface of 

a cylindrical obstacle whose axis is parallel to the magnetic field. The particle 

density and magnetic field intensity within the transition sheath on the upstream 

side are double their free stream value. The thickness of the sheath is equal to 

the cylinder radius over almost all of the upstream surface and thickens slightly 

near the joining of the upstream and downstream surfaces due to the additional 

magnetic pressure, however small, in the transition sheath. The enhanced magnetic 

pressure in the transition sheath, however slight, prevents the ions from returning 

completely to their initial free stream motion on the downstream side of the 

cylinder. As a result, a plasma stream, even though the particles initially 

have a unique velocity, will expand into the vacuum behind the cylinder. In the 

1/2 

case of a plasma with some thermal velocity (kT/2jito) the vacuum downstream 

1/2 

would be a diffuse isoceles triangle whose height would be v q (2mn/kT) R 
where v q is the initial stream velocity and R is the cylinder radius (1/2 the 
triangle base.) 

The enhanced magnetic field in the transition sheath creates a charge 
separation layer at the boundary of the undisturbed plasma stream. The total 
change in electrical potential across this layer is equal to half the stream 
pressure on the transition sheath. The plasma stream moves through the charge 
separation layer, and, on entering the transition sheath, the electrons acquire 
all the kinetic energy the ions would have had perpendicular to the magnetic field 
in the absence of the charge separation layer. At the stagnation point the 
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electron kinetic energy in the transition sheath equals the initial ion kinetic 
energy in the free stream ; but elsewhere in the transition sheath the electron 
kinetic energy equals the difference between the ion free stream kinetic energy 
and the ion stream energy in the transition layer. 

Astrophysical Applications 

The application of this work to Type I comet tails is straight forward and 
will be reported elsewhere. The slight ionization of massive cometary molecules 
through charge transfer with protons in the solar wind causes the interplanetary 
magnetic field lines to slow down drastically in overcoming the inertia of the 
massive cometary ions initially at rest. The energetic electrons in the tran- 
sition sheath which forms on the surface of the "slowed" field lines then 
rapidly and efficiently ionizes hundreds of cometary molecules per cc . The 
energetic electrons will further rapidly diffuse out along the field lines (as 
suggested by Axford in private conversion) and pull the ions out along the field 
lines so as to form the striking ray structure known as Type I comet tail rays. 

The simple model examined in this paper also yields a good qualitative 
understanding of the phenomena observed at the edge of the earth's magnetic 
field and enables predictions about the energy of the charges in the magneto- 
sheath to be made. The sphericity of the earth's magnetosphere introduces a new 
physical dimension. In the plane containing the stagnation point and the 
magnetic field lines the particle pressure both inside and outside the outer 

boundary of the transition sheath (A. J. Dessler has suggested the name magneto- 

2 

sheath) is still proportional to cos i|r. The particles are decoupled, however, from 
the magnetic field in the direction of the field, and conservation of magnetic 
flux requires that the magnetic field intensity be approximately independent 
of distanc e^tlie stagnation point or cosijr in the case of a cylinder and to increase 
away from the stagnation point in the case of a sphere. Therefore, in this plane 
the shock boundary around a circular surface would expand faster than in the 
plane at right angles to the interplanetary magnetic field. The particle density 
in this plane away from the stagnation point is less than in the plane perpen- 
dicular to the magnetic field. The particle pressure in this plane is not 
communicated along the field lines and this introduces a complication in the 
analysis depending on the isotropy of the particle velocity discussed below. 

Free hand drawings of the field lines are illustrated in Figs. 7 and 8. 
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Another aspect of the magnetosphere problem is that turbulence on the surface 
or the sharp kinks produced in magnetic field lines where the interplanetary 
field lines are joined to the lines in the magnetosheath will invalidate the 
simple two dimensional analysis presented in this work. The particle motion in 
the magnetosheath will be more nearly isotropic in three dimensions instead of 
being isotropic in two dimensions. This causes the particle density and field 
intensity to be three times the interplanetary values and the thickness of the 
sheath near the stagnation point to be more nearly one-half of the radius of 
curvature of the magnetosphere. 

Independent of these considerations the sudden increase of the magnetic 
field intensity from the interplanetary field value to the magnetosheath value 
will cause the magnetosheath electrons to have the same momentum perpendicular 
to the field lines as the protons do. Thus the particle energy and pressure 
perpendicular to the field lines will be monopolized by the electrons. The 
particle energy and pressure parallel to the field lines will be monopolized ly 
the protons, however, since charge separation effects will cause the electron 
and proton velocities to equalize. 
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FIGURE CAPTIONS 


Fig. 1 Illustration of solar wind containing a weak magnetic field perpen- 
dicular to the plane of the figure incident on a cylinder of radius R. 
The arrows represent the velocity of the field lines (particle guiding 

centers) in the free stream v , and in the transition sheath, v . The 

A ~o «*-m 

circles and loops represent partial trajectories of the particles while 
in the sheath. 

Fig. 2 Diagram of use of constant flux condition in determining the sheath 

thickness, r^ - r^, by which particles leave a volume of unit height 

after entering the volume within the length y^. 

Fig. 3 The shock structure behind a cylinder of radius R. The length of the 

2 2 1/2 

vacuum downstream from the cylinder La/R V / (v -v ) 

Fig. 4 The coordinate system used in analyzing particle motion in a magnetic 
field parallel to the z axis whose field lines travel with a velocity 
to the left of the shock boundary and with velocity v^ to the right 
of the shock boundary. 

Fig. 5 Electron and proton trajectories initially traveling with a moving 

magnetic field which is slowed down and compressed, downstream from the 
shock boundary. 

! 

Fig. 6 The geometry of the guiding center or field line motion on both sides of 
the shock boundary for oblique incidence. 

Fig. 7 Magnetic field lines indicated by dashed lines in free interplanetary 

space and in the transition sheath in the plane containing the field 

j 

direction. 

Fig. 8 Magnetic field lines, indicated by dashed lines, viewed in projection 
from the upstream side as they diverge and slide around a spherical 
obstruction bowing out most pronouncedly in the equatorial plane. 

i 





Figure 2 - Diagram of use of constant flux condition in determining the sheath thickness, r2 - rj, by which particles 

leave a volume of unit height after entering the volume within the length y 2 * 
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lines travel with a velocity v Q to the left of the shock boundary 
and with velocity v m to the right of the shock boundary. 



Figure 5 - Electron and proton trajectories initially traveling with a moving magnetic 
field which is slowed down and compressed, downstream from the shock boundary. 





Figure 6 - The geometry of the guiding center or field line motion on both sides 
of the shock boundary for oblique incidence. 





Figure 8 - Magnetic field lines, indicated by dashed lines, viewed in projection from the upstream side as they diverge 
and slide around a spherical obstruction bowing out most pronouncedly in the equatorial plane. 


